Theory of DNA translocation through narrow ion channels and nanopores with 

charged walls 
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Translocation of a single stranded DNA through genetically engineered a-hemolysin channels 
with positively charged walls is studied. It is predicted that transport properties of such channels 
are dramatically different from neutral wild type a-hemolysin channel. We assume that the wall 
charges compensate the fraction x of the bare charge qi, of the DNA piece residing in the channel. 
Our prediction are as follows (i) At small concentration of salt the blocked ion current decreases 
with x. (ii) The effective charge q of DNA piece, which is very small at x = (neutral channel) 
grows with x and at x = 1 reaches qb- (iii) The rate of DNA capture by the channel exponentially 
grows with x. Our theory is also applicable to translocation of a double stranded DNA in narrow 
solid state nanopores with positively charged walls. 



I. INTRODUCTION 

A DNA molecule in a water solution carries negative 
charge. With the help of applied voltage V, it can 
translocate through an ion channel located in a lipid 
membrane or through a solid state nanopore in a semi- 
conductor film. In this paper we are interested in the 
cases when DNA barely fits into a narrow pore leaving 
only a small gap for water with the width d < Is, where 
Ib = e 2 / HiksT is the Bjerrum length and k is the di- 
electric constant of water. An intensively studied exam- 
ple is the translocation of a single stranded DNA (ss- 
DNA) molecule through an a-hemolysin (a-HL) chan- 
nel tO, S i 11 i, 1 0, 1 1- With the average in- 
ternal diameter ~ 1.7 nm the channel can accomodate 
ssDNA molecule with ~ 1 nm diameter. In this case, 
d = 0.35 nm< Ib- 

Our theory also should be applicable to a double he- 
lix DNA (dsDNA) with ~ 2nm diameter translocating 
through a narrow solid state nanopore with 2a ~ 3 nm 
diameter [ToL [Tlj . On one hand, no experimental data is 
available for that narrow nanopores. On the other hand, 
there is impressive progres s in making and studying wider 
nanopores 0, HO, El M • 




FIG. 1: The cross-section of the membrane and the channel 
with radius a. Positive charges of the internal wall of the 
channel are shown schematically. The captured DNA is shown 
as a cylinder with radius r. A piece of DNA fitting into the 
channel is shaded. DNA phosphates of this piece have the 
total charge qb = — 12e. They are neutralized by the wall 
charges x\qb\ and the charge of mobile cations (1 — x)|gi,|. 



The peculiarity of narrow channels is related to the 
fact that the dielectric constants of the channel stem and 
lipids and the dielectric constant of the body of DNA are 
much smaller than dielectric constant of water. When 
the water filled gap between DNA and the channel wall is 
narrow (d < Ib), the electric field of small ions is squeezed 
in the gap. Potential of interaction of charges becomes 
logarithmic. This creates the electrostatic barrier for the 
ion current similar to one which was intensively studied 
for the ion transport through narrow DNA free chan- 
nels [illiEli^. 

Previous discussion of the role of this barrier for ion 
transport in the case of DNA translocation [l8[ was 
narrowly focused on neutral channels, because the wild 
type a-HL channel can be considered practically neu- 
tral. A good measure of neutrality of a channel is its 
cation/anion selectivity, measured by the ratio of cation 
to anion currents. For the wild type a-HL channel this 
ratio is equal 1.1, while it is equal to unity for a exactly 
neutral channel. It is known that at pH 7 the a-HL 
channel is neutral in the body of stem, but has the ring 
of 7 negative charges near the bottom of the the narrow 
cylindrical part of the channel (stem) [l9( . These charges 
are screened by the salt in the water outside the stem 
and, therefore, do not determine the transport through 
the channel which remains weakly selective. Solid state 
nanopores can be neutral, too. For a neutral narrow 
channel Ref. [IH addressed several challenging problems 
posed by the experimental data P, 0, 0, 0,11, [g, 0, 0, H|. 

First, Ref. [l8j explained how the electrostatic bar- 
rier makes the blocked by DNA current Ib at least 10 
times smaller than the open current Iq. Second, the ef- 
fective stall charge q s of the piece of DNA residing in 
the channel was calculated. This charge determines the 
force F s = q s V/L stalling DNA against the voltage V (L 
is the length of the channel). It was shown that for a 
neutral channel with small h/Io the charge q s ~ q^I^/I^, 
where qb is the bare charge of the piece of DNA occupying 
the channel (for a-HL channel qb = — 12e). In agreement 
with experiments this results in very small absolute value 
of the stall charge, namely q s ~ — le. Third, the origin 
of an the exponentially small and growing with the salt 
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concentration DNA capture rate was elucidated. 

Recentl y, g enetically modified a-HL channels became 
available [201 l2ll]. In this paper we concentrate on those 
of them, which have amino acids with positive residues 
on the internal wall of the narrow cylindrical part of the 
channel (stem). Internal walls of solid state nanopores 
may also be charged. The charge density of these walls 
can be tuned by different chemical treatments or just 
by a change of the solution pH. Thus, our theory for 
charged a-HL channels simultaneously addresses narrow 
charged nanopores used for dsDNA translocation exper- 
iments [loj . 

We assume below that the fraction x of the bare charge 
<7i, of DNA piece fitting into the channel is compen- 
sated by positive internal wall charges, which are roughly 
speaking randomly distributed on internal wall of the 
channel (Fig. [T]). We predict below that this simple as- 
sumption leads to a number of dramatic changes of DNA 
translocation in comparison with a neutral channel. Let 
us list these predictions: 

(i) The blocked ion current h becomes even smaller 
than in the neutral channel particularly at small concen- 
trations of salt. 

(ii) The effective charge q s of the piece of DNA resid- 
ing in the channel grows with x as q s ~ xqb- At x = 1 
the stall and bare charges of DNA are almost equal. The 
large effective charge will makes possible DNA manipu- 
lation with the help of small voltages. 

(iii) The barrier for the DNA capture decreases with 
x. As a result the DNA capture rate exponentially grows 
with x and the number of translocation events observed 
in a given experiment increases. This should lead to much 
more effective averaging of the noise and may prove help- 
ful in attempts of DNA sequencing. At some x = x c the 
capture barrier vanishes. At x c < x < 1 DNA is at- 
tracted to the channel. The capture rate then is only 
diffusion limited and independent on x. On the other 
hand, for a captured DNA the probability to escape from 
the channel becomes activated. The escape barrier grows 
with x at x > x c . 

The structure of our paper is simple. It consists of 
three sections leading to conclusions (i), (ii) and (iii) re- 
spectively. 



II. RELEASE OF COUNTERIONS AND 
BLOCKED ION CURRENT 

In the case of a-HL channel we assume the ssDNA 
molecule is a rigid cylinder coaxial with the channel. The 
inner radius of the a-HL channel is a ~ 0.85 nm, and 
the radius of the ssDNA molecule is r ~ 0.5 nm. Salt 
ions are located in the water-filled gap between them, 
with thickness 6~0.35nm. The length of the channel is 
L ~ 5 nm. This kind of model is even more appropriate 
for double helix DNA in a wider (say 4 nm in diameter) 
cylindrical solid state nanopore [IJJ, [TlJ . 

The dielectric constant of the channel or the ssDNA 



molecule (k' ~ 2) is much smaller than that of water 
(k~80). So if ssDNA is neutralized by cations and there 
is an extra charge e in the thin water-filled gap between 
the channel internal wall and ssDNA, the electric field 
lines starting from this charge are squeezed in the gap . 
This results in a high self energy of the charge [TBI, fl~6l.Tl8j] . 
According to the estimate of Ref [HI for the case of ss- 
DNA in a-HL channel, the self energy of the charge in the 
middle of the channel is ~ 5ksT. Here and everywhere 
in this paper T is the room temperature. 

Because of the large self energy of a charge in the 
narrow water gap, the piece of ssDNA inside the wild 
type neutral a-HL channel is neutralized by counterions, 
say K + in KC1 solution ssDNA covered by cations 
presents a conducting DNA backbone wire responsible for 
the blocked ion current lb at small concentration of salt 
c < 1M. In this range of concentrations lb is practically 
c independent [§]. At larger concentration c > 1M addi- 
tional pairs of anions and cations in the channel provides 
a parallel to DNA backbone wire mechanism of conduc- 
tivity. (Recall that DNA backbone wire occupies only a 
small fraction of the water filled gap.) The linear growth 
of lb with c at c > 1M is an experimental evidence for 
the second mechanism of conductivity Q 

In a mutated positively charged channel situation is 
rather different. Let us consider the channel with uni- 
formly distributed 12 positive charges (x = 1). We argue 
that in this case both ssDNA and internal wall charges 
release their counterions into the surrounding salt solu- 
tion. The net charge of the channel is still zero and, thus, 
there is practically no price in the Coulomb energy. On 
the other hand, counterion release leads to the large gain 
in their entropy. As a result DNA backbone wire looses 
its carriers and becomes an insulator. Therefore, lb is de- 
termined only by contribution of additional pairs of salt 
ions. This should lead to linear dependence of lb on c in 
the whole range of salt concentrations. In other words, 
lb becomes much smaller than in the wild type channel 
at small c < 1M, but is not strongly changed at larger 
concentration of salt. 

So far we talked about wall charges totally compensat- 
ing the bare charge of DNA (x = 1 ). At x < 1 DNA 
counterions are only partially released and conductance 
of DNA backbone wire is only partially depleted. Al- 
though the number of counterions on the DNA wire is 
proportional to x, their mobility may somewhat grow 
with decreasing x due to the increase of the number 
of empties. It is possible, but seems unlikely that this 
growth leads to the conductance maximum at x ~ 0.5. 



III. EFFECTIVE CHARGE OF DNA 

As we mentioned above for the wild type channel the 
stall charge of DNA q s is much smaller than the bare 
charge of DNA qb- Let us remind why this happens. 
Counterions neutralizing DNA in the channel receive 
from electric field momentum with the same absolute 
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value as DNA, but in the opposite direction. Most of the 
time counterions are bound to DNA charges and trans- 
fer all received momentum to DNA. During this time, 
the net electric field force acting on DNA vanishes. At 
rare moments when counterions get free and move along 
the channel contributing to If,, they transfer half of their 
momentum to the internal channel wall. This deficit of 
momentum transfer to DNA results in a small net aver- 
age force on DNA and its small effective charge q s flU ]. 

In a channel where positively charged walls compen- 
sate the bare charge of DNA, the balance of forces is 
completely different. When countrerions of DNA and 
walls are released, electric field provides opposite mo- 
mentums to DNA and to the wall charges. The latter 
are static and, therefore, transfer all their momentum to 
the wall. Thus, DNA gets its momentum only directly 
from electric field. This means that q s — qb- 

So far we talked about the channel which totally com- 
pensates charges of DNA [x = 1). Similar logic leads to 
the result q s = xqb for any x < 1. 

IV. DNA CAPTURE AND ESCAPE RATES 

Besides the blocked current and the stall charge one 
can measure the average time between the two successive 
translocation events, r, or the capture rate R c = 1/r of 
a DNA molecule into the channel. It is natural to com- 
pare the observed value of R c with the diffusion limited 
rate Rd of ssDNA capture. For the wild type neutral 
a-HL channel this comparison shows that R c <C Rd- 
The capture rate at zero voltage R c (0) is so small that 
all experiments are actually done with a large applied 
voltage V = 50 — 200 mV. Apparently there is a large 
barrier for the DNA capture. A large part of this bar- 
rier is due to the loss of the conformational entropy of 
ssDNA. The capture barrier, however, depends on the 
salt concentration, what means that a part of it has an 
electrostatic origin. The reason for such an electrostatic 
barrier is as follows [Tsi |. When a DNA molecule enters 
the channel, the DNA counterions are squeezed in the 
narrow water-filled space surrounding the DNA. Due to 



this compression the total free energy of DNA and ions is 
higher for DNA in the channel than for DNA in the bulk. 
In agreement with experiment this barrier decreases with 
growing c because the entropy of counterions in the bulk 
solution decreases and, therefore, the price for compres- 
sion is smaller. 

In the case of a channel with positively charged walls 
the ssDNA does not need to bring all its counterions 
into the channel, because there are already some posi- 
tive charges. Thus, the charge xqb is released by DNA to 
the bulk of solution making the electrostatic barrier for 
DNA smaller. Additional, roughly speaking, equal gain 
is provided by release of counterions of the wall charges, 
which screen the walls in the absence of DNA. Thus, 
due to the counterion release the electrostatic barrier be- 
comes 1 — 2x times smaller. At x = 1/2 the electrostatic 
barrier vanishes, but the conformation barrier remains 
intact and the capture rate is still activated. At x > 1/2 
the electrostatic contribution to the total barrier becomes 
negative and at small enough concentration of salt, when 
x = x c < 1, it eventually compensates the conformation 
barrier, so that R c — Rd- At x > x c the capture rate 
saturates at Rd, but the escape rate has an activation 
energy. The linear dependence of the barrier on x can be 
measured. 

To summarize, in this paper we studied DNA translo- 
cation through narrow channels with positively charged 
walls. Our predictions for the stall effective charge and 
capture rate are dramatically different from the case of 
neutral channels. Interpretation of the stall effective 
charge theory also becomes much simpler. This paper 
extends theory of narrow channels started in Ref. [18j |. 
Meanwhile, the detail, theoretical description of wider 

channels became available [Hj]. Therefore, one can 

ask when our theory crossovers to results [12| . The 
crossover happens when the width of water gap d be- 
comes larger than Ib ■ For double helix DNA this happens 
for nanopores with diameter 3.5 nm or larger. 

We are grateful to A. Aksimentiev, T. Butler, J. Gund- 
lach, Jingshan Zhang, J. J. Kasianovich, O. V. Krasil- 
nikov and A. Mcllcr for useful discussions. 
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